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Abstract 

Efficient swimming at low Reynolds numbers is a major concern 
of microbots. To compare the efficiencies of different swimmers we 
introduce the notion of "swimming drag coefficient" which allows for 
the ranking of swimmers. We find the optimal swimmer within a 
certain class of two dimensional swimmers using conformal mappings 
techniques. 

Motivation: Swimming at low Reynolds numbers is the theory of the 
locomotion of small microscopic organisms [UElEllllElinilZllHlin]- It is 

also relevant to the locomotion of small robots ^U]. Although microbots do 
not yet exist, they are part of the grand vision of nano-science [TTl IT^ IT^ 
and it is important to understand the physical constraints that underline 
their locomotion. Microbots must swim much faster than bacteria if they 
are to interface with the macroscopic world. A micron-size robot swimming 
100 times as fast as a bacterium, at the modest speed of 1 mm per second, 
has Reynolds number Re = pUL/fj, = O(10~^), and, since power scales like 
U"^, consumes 10^ more power than a bacterium. Microbots must therefore 
attempt to swim as effectively as possible, and the problem we address is 
how to search for effective swimming styles. 

Microscopic organisms use a variety of swimming techniques: Amoeba 
make large deformations of their bodies, e-coli beat flagella, paramecia use 
cilliary motion, and cyanobacteria travelling surface waves j3]. One of our 
aims is to formulate a criterium that can be used to compare different swim- 
ming styles and strokes. 

1 



We show, in the context of a two dimensional model reminiscent of 
Amoeba swimming, how one can find the optimal swimmer in a class of 
swimmers. A movie of the optimal swimmer, can be viewed in Our 
notion of optimality is closely related to a notion of efficiency which has been 
extensively used in the locomotion of microorganisms ^1 El but is more 
general and is applicable also to swimmers whose shape changes substantially 
during the swimming stroke. It is different from a notion of optimality intro- 
duced by Shapere and Wilczek P , though the two notions become equivalent 
when the amplitudes of the stroke is constrained to be small. However, as 
we shall see, small strokes are never optimal. 

The theoretical framework: Swimming results from a periodic change of 
shape. We first need to recall [S] what is meant by a shape and a located- 
shape. A located-shape is a closed surface in three dimensions (or a closed 
curve in two dimensions). The surface is parameterized so each point is 
marked and can be identified with a specific point of a fixed reference, see 
fig. El A shape is an equivalence class of congruent located-shapes that differ 
by translation and rotation. The space of all shapes consists of all such 
equivalence classes. It is an infinite dimensional space with a non-trivial 
topology. There is no a-priori metric on the space of shapes but, as we shall 
see, dissipation can be used to define a natural metric on it. 

A swimming stroke is a closed path in the space of shapes but, in general, 
an open path in the space of located-shapes. We denote the latter 7(t), 
< t < r. r is the period of the stroke. When a stroke is small the shape of 
the swimmer throughout the stroke changes only a little. Once the swimmer 
has completed a stroke, it is back to its original shape except that it is 
translated by X{'~f) and rotated. In the problem we consider the rotation 
vanishes by symmetry. 

To compute the swimming step, ^(7), and the dissipation D{'~f) associ- 
ated to the stroke 7 one needs to solve the (incompressible) Stokes equations 
for the velocity field v of the ambient liquid: 

nAv = Vp, V ■ = (1) 

subject to the boundary conditions that v vanishes at infinity and satisfies a 
no-slip condition on the surface of the swimmer. The no slip condition relates 
the liquid flow v to the swimmer movement. The latter has two parts: One 
comes from the rate of change of shape and one comes from the locomotion. 
Using internal forces, the swimmer directly controls only its shape. The 
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locomotion is determined from the requirement that at all times the total 
force and torque on the swimmer vanish 

The two dimensional case has certain special features related to the Stokes 
paradox. Specifically the condition that the total force vanishes is satisfied 
automatically and needs to be traded for the condition that a regular solu- 
tion of Eq. (PJ satisfying the boundary conditions exists, which only in two 
dimensions is not automatic. 

Optimal swimming: Optimal swimming comes from minimizing the en- 
ergy dissipated per unit swimming distance, D{'y)/X{'y), while keeping the 
average speed X{^) /r fixed. (By swimming sufficiently slowly one can always 
make the dissipation arbitrarily small.) Since both the dissipation per unit 
length, Z}(7)/X(7), and the velocity, X{pf)/T scale as 1/r a measure for the 
inefficiency of the stroke which is invariant under scaling of r is 

„ > (j(7)/-y(7)) 
'<^'> = M^wM ■ 

The smaller 5(7) the more efficient the swimmer. We call 5(7) the swimming 
drag coefficient. 5(7) is a dimensionless number in two dimensions and differs 
from the usual drag coefficient ^H] by a factor of Re. In d dimensions 5 has 
dimension L''"^. This means that (geometrically) similar swimmers, have 
the same efficiency in two dimensions, while in three dimensions smaller 
swimmers are more efficient. 

5 reduces to the notion of efficiency used in the studies of flagellar loco- 
motion pp (up to numerical and geometric factors). It is, however, different 
from a notion of optimality introduced by Shapere and Wilczek jU], where 
the dissipation per unit length is minimized while keeping the stroke period, 
rather than the velocity, flxed. 

Let I7I be the length of the stroke 7 in shape space measured using some 
metric. When a stroke is small, both ^(7) and .0(7) scale like I7P, inde- 
pendently of the choice of metric. It follows that 5(7) diverges like I/I7P for 
small strokes. Small strokes are therefore inefficient. This is in contrast with 
the Shapere- Wilczek criterion which determines the optimal stroke only up 
to an overall scale and does not penalize small strokes. 

A model swimmer with a finite dimensional shape space: Shapere and 
Wilczek jH] introduced a class of soluble models in two dimensions with a 
finite dimensional shape space. The swimmer we consider is given by the 
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image of the unit disc, \(\ < 1, under a Riemann map 
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As C traces the unit circle z = x + iy traces the boundary of the (lo- 
cated) shape in the complex plane (see Fig. [T] for examples). The space 
of located-shapes is four dimensional with coordinates {W, X,Y, Z}. X is 
naturally interpreted as the position of the swimmer, since {W, X, Y, Z} and 
{W, 0, Y, Z} describe congruent curves that differ by translation by X. Sim- 
ilarly, {W,Y,Z} and e'^'^{W,Y, Z}, G M, describe congruent closed curves 
that differ by a rotation by 0. The shape space of the model is a space of 
three complex parameters defined up to a global phase. 

When Z = and \W\ ^ \Y\ the shape is an ellipse. (When \W\ = \Y\ 
the ellipse degenerates to an interval.) A symmetry argument shows that an 
elliptic swimmer can turn but can not swim. In this sense the model with 
Z 7^ is a minimal model of a swimmer. 

For the sake of simplicity we shall, from now on, restrict {W, Y, Z} to be 
real. The shapes in this space are symmetric under mirror reflection. (This 
follows from z{(] t) — X = z{(] t) — X where z denote the complex conjugate 
of z.) A swimmer that maintains its reflection symmetry during the stroke 
can not turn and can only swim in the x direction. Hence, without loss, X(t) 
may be taken to be real as well, and the space of shapes of Eq. Q with 
real parameters, can then be identified with the three dimensional Euclidean 
space {W, F, Z} G M^. 

The solution of the model: The stroke 7 is a (parameterized) closed path 
in R^, i.e. 7 = {{W{t),Y{t), Z{t))\0 < t < r}. It generates a flow in 
the fluid surrounding the swimmer, which fills the domain corresponding to 
Id > 1. The solution of the Stokes equations, Eq. ((H), can then be obtained 
by conformal mapping methods [3[ '■ 
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where dot denotes time derivative. The flow vanishes at inflnity provided 
/2(oo) = /i(oo) = 0. From Eq. (jH)) one flnds that this is the case provided: 

dX = AdY, A=— £— . (7) 

This is the basic relation between the swimming (the response), dX, and the 
change in shape (the controls), {W, Y, Z}. The notation dX stresses that X 
does not integrate to a function of {W, Y, Z}. Geometrically, this relation is 
interpreted as a connection on the space of shapes Note that A{W, Y, Z) 
is a homogeneous function of degree zero. 

The power P dissipated by the swimmer is calculated by integrating the 
stress times the velocity on the surface of the swimmer: 

P = Im ^ V ^fi ^^^^^ dz — pdz^ , (8) 

where p is the fluid pressure, p = — 4/iRe/{(2). Using the explicit solution 
given by Eqs. (jUSEl), one obtains 

P = 47r/i (W^2 + y2 + _ 
The dissipation of a stroke is then 

D(7)=47r/i/" {W^ + Y^ + Z^)dt . (10) 
Jo 

The Physical cone: A physical shape does not self-intersect. Since there 
are points {W, Y, Z} G which represent curves that self-intersect, e.g. Z = 
0, W = ±y, we need to remove them. The physical shapes make a cone, for 
if {W, Y, Z} does not self-intersect neither does X{W, Y, Z} with A > 0. The 
shapes associated with points in the interior of the cone are smooth. Points 
on its boundary correspond to shapes with cusp-like singularities (precursors 
of self-intersections). The boundary of the physical cone is given by those 
{W, Y, Z} for which there exists ( of unit modulus such that z'{() = 0. One 
flnds that the physical cone around the W axis is bounded by two planes and 
a quadratic surface and is given by , see Fig. [TJ 

W>Y±V2Z, WY>2Z^-W^. (11) 

Optimization and orbits in magnetic fields: Admissible strokes are closed 
paths 7 that lie in the physical cone. Consider the problem of minimizing 
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Figure 1: Each point {W,Y,Z} in the cone on the left corresponds to a two 
dimensional curve that does not self intersect shown on the right. The colors 
of the shapes (right) matches the color of the dots (left). The two ellipses 
correspond to the two interior points. The six shapes with cusps correspond 
to the six points on the boundary of the physical cone. 



the dissipation, -D(7), of Eq. Q subject to the constraint that the step size 
is 

X(7) = (f AdY. (12) 

This is a standard problem in variational calculus. Note that since one may 
set r = 1 without loss, fixing the average speed is equivalent to fixing the 
step size X{'j). The minimizer, 7(t), must then either follow the boundary, 
or solve the Euler-Lagrange equation of the functional 

Sq{-i) = A'Kii I {W^ + Y'^ + Z^)dt + q AYdt, (13) 
Jo Jo 

where g is a Lagrange multiplier. Sg can be interpreted as the action of a 
non-relativistic particle whose mass is Svr/i and whose charge is q moving in 
three dimensions under the action of a magnetic field with vector potential 
AY. X(7) may be interpreted as the fiux of B{W, Y, Z) = V x (AY) through 
the closed path 7. 

The parameterized stroke 7(t) that minimizes Sg has constant velocity 
I7I = const. (This follows from the fact that the fiux ^(7) is independent of 
parametrization and that the action of a free particle along a one dimensional 
curve is minimized at constant speed.) The dissipation is then .0(7) = 
47r/i— , where I7I is the length of the orbit, and the drag 5(7) is simply: 
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Since the variational problem is for a domain with non-smooth bound- 
aries, Fig. Q one may worry if the minimizer fails to be smooth. This is not 
the case. For if 7(t) has a corner, smoothing the corner on a scale e shortens 
the length I7I by 0{e) while |X(7)| varies by 0{e^). In particular, it fol- 
lows that the minimizer avoids the corners of the physical cone. Moreover, 
whenever it hits or leaves the boundary of the physical cone, it does that 
tangent ially, without corners. 

Incompressible swimmers: Incompressible swimmers make a natural (and 
biologically important) class of swimmers. In two dimensions incompressibil- 
ity implies constant area. The area of the swimmer whose shape is given by 
Eq. © is 



Fixing the area of the swimmer corresponds to restricting the stroke to a 
hyperboloid in shape space. We choose the unit of area so that the area of 
the swimmer is vr. The intersection of the constant area hyperboloid with 
the cone of physical shapes is the Big-Ben shaped surface shown in Fig. |21 
Physical strokes are represented by closed paths that lie inside this domain, 
and our aim is to find the stroke that minimizes the swimming drag ^(7). 

Large strokes are inefficient: The model admits strokes that extend to 
arbitrarily large values of Y and W . Since, by Eqs. (|7I11I12|) . the total flux 
of B through the Big-Ben shaped region of Fig.|21is infinite, the swimmer can 
swim arbitrarily large distances with a single stroke. However, as we presently 
show, large strokes are inefficient. The domain of physical incompressible 
shapes for a swimmer of area vr is contained in the strip \Z\ < 1. Since 
A = 0{y), for Y large, a long excursion, of order i, in the Y direction 
contributes 0(log£) to ^(7), but 0{i) to |7|. Therefore, as ^ ^ 00 the drag 
coefficient 6 diverges like i/\ogi. 

The minimizer: Since the drag diverges for small strokes and also for large 
strokes (for an incompressible physical swimmer) the minimizer of 5 is a finite 
stroke. It can be computed numerically using the following procedure: Since 
the minimizer is independent of the period r, one may, without loss, restrict 
oneself to orbits with fixed energy, say E = 1. Pick the charge q and find 
a smooth and closed orbit on the hyperboloid of constant area with (four) 
sections in the interior of the cone of physical shapes and (four) sections on 
its boundary. There is a unique such orbit 7^ for all q that are small enough. 
For each such orbit one computes, numerically, ^(7^), and \'jq\ to get ^(7^) 




(15) 
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Figure 2: The Big-Ben shaped surface is the intersection of the physical 
cone (no self-intersection) with the hyperboloid of shapes with fixed area. 
The boundary of the domain of incompressible non-self-intersecting shapes 
is drawn black. The inscribed red curve is the optimal stroke as computed 
numerically. 

from Eq. (|14p. What remains is a minimization problem in one variable, q, 
which yields the optimal stroke. The optimal stroke in shape space is shown 
in Fig.|21while snapshots of the corresponding swimming motion in real space 
are shown in Fig. El For the optimal stroke we find (^optimal ~ 9.12 . For the 
sake of comparison with a squirming circle consider 7 which is a small circle 
of radius r in the Y — Z plane, centered at Y = Z = 0,W = Wq. One readily 
finds from Eqs. ()7ll4j) that 6 = 8{Wo/r)'^. A squirmer that changes it shape 
by 10% has, r/Wo = 0.1, and 5 = 800. 

Perspective: The optimal swimmer we have found is optimal within the 
class of Riemann maps of Eq. Q satisfying incompressibility. Enlarging the 
class of Riemann maps, would allow for better swimmers. It is conceivable 
that there are superior swimmers that use quite different swimming styles. 
The importance of the model lies in that it demonstrates a scheme for a 
systematic search of efficient swimmers, and provides benchmark for 6 for 
better swimmers to beat. 

Acknowledgment: We thank E. Braun , G. Kosa, and D. Weihs for useful 
discussions, E. Yariv for pointing out ref. ^U] and U. Sivan for proposing the 
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Figure 3: Snapshots of the optimal swimmer shifted vertically for visibility. 
The top and bottom snapshots are related by a (horizontal) translation. The 
shapes with cusps correspond to those parts of the stroke that lie on the 
boundary of the domain of simple shapes. The (red) dots are fixed reference 
points in the body. When the swimmer is approximately triangular, the 
base of the triangle functions as an anchor that pushes or pulls the opposite 
vertex. 
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